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Abstract 



We present a study of events with Z bosons and hadronic jets produced in 
pp collisions at a center-of-mass energy of 1.8 TeV. The data consist of 6708 
Z — ► e + e~ decays from 106 pb~ 4 of integrated luminosity collected using the 
CDF detector at the Tevatron Collider. The Z + > n jet cross sections and jet 
production properties have been measured for n = 1 to 4. The data compare 
well to predictions of leading order QCD matrix element calculations with 
added gluon radiation and simulated parton fragmentation. 

PACS numbers: 14.70 
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The Z boson was discovered in pp collisions at CERN in 1983 Since then, the 

hadronic production properties of the Z have been studied using relatively small event 
samples, most recently from experiments at the Fermilab Tevatron Collider ||. In this 
Letter we describe an extension of these studies based on a much larger sample of 6708 
Z — ► e + e _ decays obtained from 106 pb -1 of integrated pp luminosity using the CDF 
detector. This event sample is large enough to study the production properties of high 
energy hadronic jets associated with Z boson production. The pp —>■ Z + jet events provide 
a good test of Quantum Chromodynamics (QCD) calculations since the event sample has 
very small background and the presence of the Z selects high Q 2 parton-level processes. A 
determination of the reliability of QCD calculations for heavy boson production is important 
for verifying the direct W boson backgrounds to top quark production Z boson events 

are free of top quark contamination and provide a clean test of Standard Model heavy boson 
production. As part of this study we examine the Z boson sample for evidence of excess b 
quark decays that could indicate new particle production. 

The elements of the CDF detector of primary importance to this analysis are the cen- 
tral tracking chamber (CTC), the calorimeters, and the silicon vertex detector (SVX). The 
CTC, which is immersed in a 1.4 T solenoidal magnetic field, measures the momenta and 
trajectories of charged particles in the region \rj\ < 1.1 (where rj = — ln(tan(6 l /2)) ||). 
The calorimeters are divided into electromagnetic and hadronic components and cover the 
pseudorapidity range \r]\ < 4.2. The four-layer silicon strip detector (SVX) @, located just 
outside the beampipe, provides precise tracking in the plane transverse to the beam direction 
and is used to reconstruct secondary decay vertices from B hadrons. The CDF detector is 
described in detail elsewhere J|,|8|]. 

Events with Z bosons are identified by looking for the decay Z — > e + e~ . We demand that 
candidate events pass a high-i?^ electron trigger. The event selection requires an isolated 
electron |J that satisfies tight selection cuts |L(| and has a track in the central region 
(|^| < 1.1). Events containing a central electron with transverse energy Ex > 20 GeV are 
required to have a second electron that satisfies looser selection cuts and has opposite charge 
where the sign of the track curvature is well-measured. The second electron is detected in 
either the central calorimeter {Ex > 20 GeV and \rj\ < 1.1), the plug calorimeter {Ex > 15 
GeV and 1.1 < \r}\ < 2.4), or the forward calorimeter {Ex > 10 GeV and 2.4 < \rj\ < 3.7). 
We remove electrons with photon conversion characteristics. The separation in t] — <p between 
the centroid of an electron calorimeter cluster and that of any jet in the event, measured in 
AR = VA?+A^, must exceed AR = 0.52. 

An event sample of 6708 Z — > e + e~ decays is selected by requiring the electron pair 
masses to be within 15 GeV/c 2 of the nominal Z boson mass of 91 GeV/c 2 . The background 
in these data and the Z — > e + e~ acceptance due to our selection cuts are measured as a 
function of associated jet multiplicity as described below. 

Hadronic jets produced in association with the Z bosons are selected using a clustering 



algorithm |IT| with a cone size of Rj = 0.4. The jet energies are corrected to account 
for variations in calorimeter response, estimates of fragmentation energy outside the jet 
cone, and underlying event energy within the jet cone. Initially, we consider all jets with 
corrected transverse energy Ex > 12 GeV. When the separation between two jets is less than 
AR = 0.52, they are combined vectorially into a single jet. The requirements Ex > 15 GeV 
and \t}\ < 2.4 are applied to obtain the hadronic jets used for this analysis. In our sample 
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of 6708 Z events, 1310 have >1 jet, 279 have >2 jets, 57 have >3 jets, and 11 have >4 jets. 
We correct these jet multiplicities for two additional effects: photons counted as hadronic 
jets, and jets produced in extra pp interactions that occur in the same bunch crossing as 
the Z event. A Z + photon Monte Carlo calculation ]12| with detector simulation yields 
a correction varying from -2% to -3% for photons as a function of jet multiplicity. Using 
minimum bias events, we estimate the number of jets from additional interactions that pass 
the selection cuts, and obtain corrections varying from -3% to -5%. 

The backgrounds to the Z boson are dominated by jets faking electrons, but include some 
contributions from heavy quark, W — > ev + jet, and Z — ► t + t~ decays. To measure the 
background we employ a data sample in which all Z boson selection cuts have been applied, 
except the mass window and electron isolation cuts. By selecting events from this sample in 
which neither electron candidate is isolated, we obtain a set of events that is almost entirely 
background with no measurable contribution from Z — > e + e~. The mass distribution for 
these background events is independent of electron candidate isolation. This allows us to 
estimate the number of background events with isolated electrons that lie within the mass 
window. This method yields background estimates in the Z boson event sample that are 
small even at high jet multiplicities. The la upper limits for the backgrounds are 1.1%, 
2.3%, 3.0%, and 4.0% for the n > 0, 1, 2, and 3 jet events, respectively. 

The acceptance for the electrons from a Z boson decay to pass all selection cuts has 
been measured as a function of the number of associated hadronic jets. Losses due to the 
electron Et cuts and detector geometric acceptance requirements have been studied using 



an inclusive Z boson production model and a leading-order Z + jets QCD calculation [13 



The loss of electrons due to overlap with hadronic jets is determined directly from the data 
by taking Z — > e + e~ events and using a Monte Carlo program to re-decay the Z bosons. 
The efficiency for decay electrons to be separated from jets and to be isolated sufficiently 
from other energy in the calorimeters is determined as a function of jet multiplicity. The 
systematic uncertainties in the overlap efficiencies are determined by varying the polarization 
of the Z and the separation between the electrons and any jets in the event. For electrons 
passing the E T , geometric acceptance, and electron-jet separation cuts, a final acceptance 
correction is made for electron identification cuts and the efficiency of the online trigger, 
which selects a high-i?y central electron. The total Z — > e + e~ detection efficiencies vary 
from (37.3 ± 0.7)% for Z bosons without jets to (31.5 ± 1.4)% for Z bosons with > 4 jets. 

The systematic uncertainties on the number of jets due to the jet selection cuts are 
determined by varying the jet energy scale by ±5%, the \rj\ cut by ±0.2, the underlying 
event correction by ±50%, the photon removal correction by ±15%, and the probability 
of jets from additional pp interactions by ±100%. The total jet-counting uncertainties are 
found to range between 11% for the > 1 jet to 23% for the > 4 jet sample. These errors 
dominate the uncertainties in the cross section measurements, except for Z ± > 4 jet events 
where the statistical error is larger. 

We measure the cross section for Z production as a function of jet multiplicity using 
the ratio of the number of detected Z events with > n jets to the total number of detected 
Z events. The cross sections are then determined from this ratio using the inclusive Z 
boson cross section of 231 ± 12 pb [T3J] and the ratios of detection efficiencies. This method 
takes advantage of the cancellation of some systematic uncertainties in the ratios, and gives 
the most accurate relative Z + > n jet cross sections. The measured cross sections for Z 
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+ > n (n = 1 to 4) jet events are given in Table pj as are the measured ratios of a(Z+ > n 
jets)/cr(Z+ > n — 1 jets), which are constant at ~0.2 for n = 1 to 4. 

The measured cross sections and the production characteristics of the hadronic jets can 
be compared to QCD calculations. The leading-order matrix element calculations for Z + n 
parton events (for n — 1, 2, 3) are obtained with the VECBOS ||13[ Monte Carlo program 
| 15| . We use both MRSA and CTEQ3M parton distribution functions, the two- loop a s 
evolution, and factorization and renormalization scales varying from Q 2 = M\ + p\ of 
the Z boson to Q 2 = (pt) 2 , where (pt) is the average px of the generated partons. The 
QCD predictions are indistinguishable for the two parton distribution functions within the 
statistical uncertainties of our calculation. 

The QCD-predicted Z + jet events are obtained from the parton-level events by including 
gluon radiation and hadronic fragmentation using the HERWIG [ |I~6| , |T7| shower simulation 



algorithm JT8|. This procedure represents a partial higher-order correction to tree- level 
diagrams. The Z boson events with hadron showers are then introduced into a full CDF 
detector simulation, and the resulting jets are identified and selected in a manner identical to 
the data. This procedure allows us to make direct comparisons between the QCD predictions 
and data. 

The Monte Carlo events passing our jet selection cuts are used to determine QCD- 
predicted cross sections that are compared to the experimental measurements corrected for 
Z — > e + e~ decay losses and backgrounds. These calculated QCD cross sections are given 
in Table |, and compared to the data via ratios of the measured to QCD-predicted cross 
sections. Figure [l] shows a plot of the measured cross sections for Z + > n jets as a function 
of n, with the QCD predictions resulting from variations in Q 2 indicated by a superimposed 
band. For Q 2 = (pr) 2 , the measured Z + > n jet cross sections range from 0.83 to 1.29 
times the leading-order QCD predictions. For Q 2 = M§ + p Tz , the ratios of cross sections 
to leading-order QCD predictions for each n are larger, but nearly constant at 1.7. 

In comparing the jet production properties of Z + jets events, we use the same Monte 
Carlo event generation described above but apply both Z — > e + e _ and jet selection cuts. 
For this study the QCD-predicted event distributions with Q 2 = (pr) 2 are normalized to 
the number of events in the data samples. Corrections are made for photons counted as 
hadronic jets and for jets from extra pp interactions in the same crossing as the Z event, but 
no corrections are made for the very small Z boson backgrounds. Figure shows the Et 
spectra for the first, second, and third jets (ordered by decreasing Et) in the > 1, > 2, and 
> 3 jet Z events, respectively. In the Z + > 1 jet event sample we measure the angle 0* 
between the Z boson and the average beam direction in the Z + leading jet center-of-mass 
frame. The distribution of the quantity | cosG*|, for the range | cosG*| < 0.95, is shown in 
Figure |3|. In events with two or more jets, we measure the separation in 7] — <p space between 
the two leading jets. The resulting ARjj distribution is given in Figure U This distribution 
is not corrected for the small number of photons counted as hadronic jets. 

The measured jet production properties agree well with the QCD predictions. The \ 2 
per degree of freedom is 16.1/18 for the cosG* distribution shown in Figure |3| and 7.8/7 
for the ARjj distribution shown in Figure £|. The QCD-predicted ARjj distribution is 
sensitive to the manner in which gluon radiation is added by the HERWIG algorithm to the 
VECBOS matrix element, and limiting this radiation results in poorer agreement with the 
data. The Et distributions shown in Figure are also in reasonable agreement with the 
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QCD predictions. The values of the \ 2 P er degree of freedom are 19.8/12, 3.7/4, and 6.2/2 
for the jet E T spectra shown in Figures ||a, ||b, and ||c, respectively. 

The Z + jets data sample has been examined for secondary vertices that are characteristic 
of b quark decays. The secondary vertices are measured using the SVX detector and are 
reconstructed using the same secondary vertex algorithm developed by CDF and used in the 
top quark search @||. Of the 1665 jets in the sample, 442 are candidates for secondary vertex 
reconstruction, which means that they have at least two tracks in the SVX, an uncorrected 
jet Et > 15 GeV, and < 2. A control sample of inclusive jet events is used to determine 
the number of secondary vertices expected from jets in QCD events with no heavy boson. 
In all, six secondary vertex candidates are found in the data sample where 6.3 ± 1.0 are 
expected, so there is no excess of b quark decays over the expected number in this Z + jets 
dataset. 

In summary, this Letter contains an analysis of jet production properties associated with 
Z — > e + e~ events selected from 106 pb _1 of pp collisions at a center-of-mass energy of 
1.8 TeV. Comparisons are made between the data and leading order parton matrix elements 
with HERWIG-simulated parton shower and fragmentation. The ratios of the measured to 
QCD-predicted cross sections are found to vary from 1.29 ± 0.17 to 1.69 ± 0.22 for pp — ► 
Z+ > 1 jet when the Q 2 scale is varied from (px) 2 to M§ + p\ . The QCD-predicted 
jet production properties are in generally good agreement with the measured distributions. 
The incidence of b quark decays in the jets associated with Z bosons is consistent with that 
observed in similar jets from events with no heavy boson. 
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Fisica Nucleare; the Ministry of Education, Science and Culture of Japan; the Natural 
Sciences and Engineering Research Council of Canada; the National Science Council of the 
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FIGURES 




FIG. 1. Cross section for Z + > n jets versus n, where Z — > e + e~. The QCD prediction band 
is determined by using two different renormalization scales. The solid line is an exponential fit to 
the measured cross sections. The error bars include statistical and systematic uncertainties. 
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FIG. 2. Transverse energy of the (a) first, (b) second, and (c) third highest Er jets in >1, >2, 
and >3 jet events, respectively. The points are the data with statistical errors only, and the solid 
histograms are the QCD predictions (normalized to the data) described in the text. 
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FIG. 3. The | cos 0*| distribution of the Z in >1 jet events. The angle 0* is the angle of the Z 
measured with respect to the average beam direction in the Z + leading jet center-of-mass frame. 
The Monte Carlo predictions are normalized to the data, and the data errors are statistical only. 
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FIG. 4. Separation in 77 — <f> space between the two leading jets in Z + > 2 jet events. The 
Monte Carlo predictions are normalized to the data. The data errors are statistical only. 
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TABLES 



n 
Jets 


BR X CTData 

(Pb) 


Q 2 = (pt) 2 


Q 2 = M 2 +p 2 Tz 


a{n)/a{n — 1) 
(Data) 


BR ■ o-QCD(pb) 


0"Data/cQCD 


BR ■ o-QCD(pb) 


CData/^QCD 


> o 


231±6±11 






223 ± 10 


1.04 ± 0.07 




> 1 


45.2 ± 1.2 ± 5.7 


35.16 ± 0.54 


1.29 ± 0.17 


26.82 ± 0.40 


1.69 ± 0.22 


0.196 ± 0.007 


> 2 


9.7 ± 0.6 ± 1.8 


10.53 ± 0.38 


0.92 ± 0.19 


5.77 ± 0.16 


1.68 ± 0.34 


0.215 ± 0.014 


> 3 


2.03 ± 0.28 ± 0.49 


2.44 ± 0.17 


0.83 ± 0.24 


1.23 ± 0.08 


1.66 ± 0.47 


0.210 ± 0.027 


> 4 


0.43 ± 0.13 ± 0.11 










0.211 ± 0.059 



TABLE I. Z + > n jet Cross Sections. The first error on the data cross sections is the statistical 
error; the second includes the systematic error on the Z acceptance and the luminosity error added 
in quadrature. For n = 1 to 4, the second error also includes the jet-counting uncertainty, as 
described in the text. The leading order (LO) Monte Carlo cross sections are generated with 
VECBOS for Q 2 scales of (pt) 2 and M§ +Pt z - The > jet calculation is next-to-next-to-leading 
order (NNLO) with Q 2 = M|. The error on the NNLO QCD cross section includes systematic 
and statistical uncertainties; for the LO cross sections the errors are statistical. 
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